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i— ■ ' The combined HERA-I data set, of neutral and charged current inclusive cross-sections 

, for e + p and e~p scattering, is used as the sole input for a next-to-leading order (NLO) 

QCD parton distribution function (PDF) fit. The consistent treatment of systematic 
CO ' uncertainties in the joint data set ensures that experimental uncertainties on the PDFs 

can be calculated without need for an increased \ 2 tolerance. This results in PDFs 
with greatly reduced experimental uncertainties compared to the separate analyses 
of the ZEUS and HI experiments. Model uncertainties, including those arising from 
parametrization dependence, are also carefully considered. The resulting HERAPDFs 
I ' have impressive precision compared to the global fits. 
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The kinematics of lepton hadron scattering is described in terms of the variables Q 2 , the 
invariant mass of the exchanged vector boson, Bjorken x, the fraction of the momentum of 
the incoming nucleon taken by the struck quark (in the quark-parton model), and y which 
■ measures the energy transfer between the lepton and hadron systems. The differential cross- 

section for the neutral current (NC) process is given in terms of the structure functions by 
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where Y± = 1 ± (1 — y) 2 . The structure functions F% and xF% are directly related to quark 
distributions, and their Q 2 dependence, or scaling violation, is predicted by perturbative 
QCD. For low x, x < 10 -2 , F2 is sea quark dominated, but its Q 2 evolution is controlled 
by the gluon contribution, such that HERA data provide crucial information on low-x sea- 
quark and gluon distributions. At high Q 2 , the structure function XF3 becomes increasingly 
important, and gives information on valence quark distributions. The charged current (CC) 
interactions also enable us to separate the flavour of the valence distributions at high-21, 
since their (LO) cross-sections are given by, 
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Parton Density Function (PDF) determinations are usually obtained from global NLO 
QCD fits [21 Ell!], which use fixed target DIS data as well as HERA data. In such analyses, 
the high statistics HERA NC e + p data have determined the low- a; sea and gluon distribu- 
tions, whereas the fixed target data have determined the valence distributions. Now that 
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high-Q 2 HERA data on NC and CC e + p and e~p inclusive double differential cross-sections 
are available, PDF fits can be made to HERA data alone, since the HERA high Q 2 cross- 
section data can be used to determine the valence distributions. This has the advantage 
that it eliminates the need for heavy target corrections, which must be applied to the ^-Fe 
and fiD fixed target data. Furthermore there is no need to assume isospin symmetry, i.e. 
that d in the proton is the same as u in the neutron, since the d distribution can be obtained 
directly from CC e + p data. 

The Hf and ZEUS collaborations have both used their data to make PDF fits [5], [§]■ 
Both of these data sets have very small statistical uncertainties, so that the contribution of 
systematic uncertainties becomes dominant and consideration of point to point correlations 
between systematic uncertainties is essential. The ZEUS analysis takes account of corre- 
lated experimental systematic errors by the Offset Method, whereas HI uses the Hessian 
method [TJ. Whereas the resulting ZEUS and HI PDFs are compatible, the gluon PDFs do 
have rather different shapes, see Fig.[TJ and the uncertainty bands spanned by these analyses 
are comparable to those of the global fits. 

It is possible to improve on this situation since ZEUS and HI are measuring the same 
physics in the same kinematic region. These data have been combined them using a 'theory- 
free' Hessian fit in which the only assumption is that there is a true value of the cross- 
section, for each process, at each x, Q 2 point [8], [9]. The resulting systematic uncertainties 
on each of the combined data points arc significantly smaller than the statistical errors. In 
the present paper this combined data set is used as the input to a NLO QCD PDF fit. 
The consistency of the input data set and its small systematic uncertainties enable us to 
calculate the experimental uncertainties on the PDFs using the \ 2 tolerance, A% 2 = 1. 
This represents a further advantage compared to those global fit analyses where increased 
tolerances of A\ 2 — 50 — 100 are used to account for data inconsistencies. 

For the present HERAPDF0.1 fit, the role of correlated systematic uncertainties is no 
longer crucial since these uncertainties are relatively small. This ensures that similar results 
are obtained using either Offset or Hessian methods, or by simply combining statistical 
and systematic uncertainties in quadrature. For our central fit we have chosen to combine 
the 43 systematic uncertainties which result from the separate ZEUS and HI data sets in 
quadrature, and to Offset the 4 sources of uncertainty which result from the combination 
procedure. This results in the most conservative uncertainty estimates on the resulting 
PDFs. 

Despite our conservative procedure the experimental uncertainties on the resulting PDFs 
are impressively small and a thorough consideration of further uncertainties due to model 
assumptions is necessary. In section [2] we describe the NLO QCD analysis and model 
assumptions. In section |3] we give results and in section [4] we give a summary. 

2 Analysis 

The QCD predictions for the structure functions are obtained by solving the DGLAP evo- 
lution equations at NLO in the MSbar scheme with the renormalisation and factorization 
scales chosen to be Q 2 . The DGLAP equations yield the PDFs at all values of Q 2 provided 
they are input as functions of x at some input scale Qq. This scale has been chosen to be 
Qq = 4GeV 2 and variation of this choice is considered as one of the model uncertainties. 
The resulting PDFs are then convoluted with NLO coefficient functions to give the structure 
functions which enter into the expressions for the cross-sections. The choice of the heavy 
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quark masses is, m c = 1.4, m& = 4.75GeV, and variation of these choices is included in the 
model uncertainties. For this preliminary analysis, the heavy quark coefficient functions 
have been caluclated in the zero-mass variable flavour number scheme. The strong coupling 
constant was fixed to a s (Mz) = 0.1176 [10], and variations in this value of ±0.002 have also 
been considered. 

The fit is made at leading twist. The HERA data have of the invariant mass of the 
hadronic system,^ 2 , of W^ in = 300GeV 2 and maximum x, x max = 0.65, such that they 
are in a kinematic region where there is no sensitivity to target mass and large-x higher 
twist contributions. However a minimum Q 2 cut is imposed to remain in the kinematic 
region where perturbative QCD should be applicable. This has been chosen such that 
Qmin — 3-5 GeV 2 . Variation of this cut is included as one of the model uncertainties. 

A further model uncertainty is the choice of the initial parametrization at Qq. The PDFs 
are parametrized by the generic form 

xf(x) = Ax B (l - x) c (l + Dx + Ex 2 + Fx 3 ), (1) 

and the number of parameters is chosen by 'saturation of the x 2 \ such that parameters 
D, E, F are only varied if this brings significant improvement to the \ 2 ■ Otherwise they are 
set to zero. 

For our central fit, the PDFs which arc parametrized arc xu v , xd v , xg and xU, xD. 
The normalisation parameters, A, for the d and u valence are constrained to impose the 
number sum-rules and the normalisation parameter A for the gluon is constrained to impose 
the momentum sum-rule. The B parameters which constrain the low-x behaviour of the 
u and d valence distributions are set equal, and the B parameters are also set equal for 
xU and xD, such that there is a single B parameter for the valence and another different 
single B parameter for the sea distributions. Assuming that the strange and charm quark 
distributions can be expressed as x independent fractions, f s = 0.33 and f c = 0.15, of 
the d and u type sea, gives the further constraint A(U) = A(D)(l — / s )/(l — / c ). The 
value of f s = 0.33 has been chosen to be consistent with determinations of this fraction 
using neutrino induced di-muon production. This value has been varied to evaluate model 
uncertainties. The charm fraction has been set to be consistent with dynamic generation of 
charm from the start point of Q 2 = m 2 in a zero-mass-variablc-fiavour-number scheme. A 
small variation of the value of f c is included in the model uncertainties. Saturation of the 
X 2 leads us to set the parameters D, E, F = 0, for all partons except xu v for which only 
F = 0. 

The results are presented using this parametrization, including six sources of model 
uncertainty due to variation of: m c , mi,, f s , f c ,Qo, Q 2 m in- Comparison is made to three other 
classes of parametrization, one based on the ZEUS-JETs parametrization [5], one based 
on the HI parametrization [6] and one based on the current parametrization but allowing 
O / for the gluon. Comparison is also made to results obtained by varying a s (Mz), see 
reference [T] for details. Our central choice has less model dependence than the ZEUS-Style 
parametrization because it has fewer asumptions concerning d — u, and it has less model 
dependence than the Hl-style parametrization in that it does not assume equality of all B 
parameters. Furthermore, although all types of parametrization give acceptable \ 2 values, 
the central parametrization has the best \ 2 an( i it gives the most conservative experimental 
errors. 
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HI and ZEUS Combined PDF Fit 
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Figure 1: Left: PDFs from the ZEUS- JETS and H1PDF2000 PDF separate analyses of 
ZEUS and HI. Right: HERAPDF0.1 PDFs from the analysis of the combined data set 

3 Results 

The total uncertainty of the PDFs obtained from the HERA combined data set is much 
reduced compared to the PDFs extracted from the analyses of the separate HI and ZEUS 
data sets, as can be seen from the summary plots in Fig. [IJ where these new HERAPDF0.1 
PDFs are compared to the ZEUS-JETS and H1PDF2000 PDFs. 

In Fig.[2]we show the HERAPDF0.1 PDFs compared to the CTEQ6.1 PDFs, which also 
use a zero-mass variable flavour number scheme, and to the preliminary MSTW08 PDFs [11] , 
which use a massive variable flavour number scheme. The precision of the HERAPDF0.1 
for the low- a; sea and gluon is impressive. 



4 Summary 

Now that high-Q 2 HERA data on NC and CC e + p and e~p inclusive double differential cross- 
sections are available, PDF fits can be made to HERA data alone, since the HERA high 
Q 2 cross-section data can be used to determine the valence distributions and HERA low Q 2 
cross-section data can be used to determine the Sea and gluon distributions. The combined 
HERA-I data set, of neutral and charged current inclusive cross-sections for e + p and e~p 
scattering, has been used as the sole input for a NLO QCD PDF fit in the DGLAP formal- 
ism. The consistent treatment of systematic uncertainties in the joint data set ensures that 
experimental uncertainties on the PDFs can be calculated without need for an increased x 2 
tolerance. This results in PDFs with greatly reduced experimental uncertainties compared 
to the separate analyses of the ZEUS and HI experiments. Model uncertainties, including 
those arising from parametrization dependence, have also been carefully considered. The 
resulting HERAPDFs have impressive precision compared to the global fits. 
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